Evolution of spectral properties along the 0(6)-U(5) transition 
in the interacting boson model. II. Classical trajectories 
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This article continues our previous study of level dynamics in the [0(6)— U(5)]dO(5) transition of 
the interacting boson model nucl-th/0504016 using the semiclassical theory of spectral fluctuations. 
We find classical monodromy, related to a singular bundle of orbits with infinite period at energy 
E = 0, and bifurcations of numerous periodic orbits for E > 0. The spectrum of allowed ratios of 
periods associated with 0- and 7-vibrations exhibits an abrupt change around zero energy. These 
findings explain anomalous bunching of quantum states in the E « region, which is responsible 
for the redistribution of levels between 0(6) and U(5) multiplets. 
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I. INTRODUCTION 

In the first part of this work p| (hereafter referred to 
as Part I), we have discussed the evolution of level en- 
ergies and wave functions along the [0(6)— U(5)]dO(5) 
transition in the interacting boson model (IBM) |2|. It 
is known that this transitional class is integrable — due 
to the 0(5) underlying symmetry and the associated se- 
niority quantum number v — and exhibits a second-order 
ground-state phase transition from deformed 7-soft to 
spherical equilibrium shapes. 

Remind that our family of model Hamiltonians is given 

by 



H(rj) = a 



(i) 



with 77 € [0,1] denoting a dimensionless control param- 
eter that drives the system between the 0(6) (77 = 0) 
and U(5) (77 = 1) dynamical symmetries. The spectrum 
of H(rj) at any point of the transitional path depends 
on the specific interplay of both terms in Eq. , where 
Q = [«t d+ eft s]( 2 ) represents the quadrupole operator and 
hd = ((P ■ d) the d-boson number operator. Note that 
N is the total number of bosons, which in the classical 
limit tends to infinity (both terms in the above Hamilto- 
nian are properly normalized by the N k denominators to 
yield finite contributions in this limit), and a — 1 MeV 
is an arbitrary scaling factor. The N — > 00 ground-state 
shape-phase transitions takes place at rj c — |. 

We have shown that one of the most significant features 
of spectra in the 77 g [0, 1] transitional regime of Hamil- 
tonian is the pattern of alternating compressions and 
dilutions of levels with angular momentum I — around 
energy E w 0. This pattern spreads over a wide interval 
of the control parameter between 77 w 0.3 and 0.8, see 
Fig. 1 in Part I. After deconvoluting spectra with differ- 
ent seniorities, it transforms into a sequence of avoided 
crossings that constitute what we called the "shock-wave 
scenario" 0. 



The level bunching pattern represents basically a huge 
oscillation of the level density in the E w region, not 
dissimilar to shell effects in single-particle spectra of some 
quantum- mechanical potentials. There exists a deep and 
far-reaching relation between fluctuations of the quan- 
tum level density and properties of periodic orbits in the 
classical counterpart of the given system [j, 3 ■ While it 
is known that each periodic orbit brings one oscillatory 
term into the level density, with an amplitude related to 
the orbit's dimensionality and stability [a, |6J, |j| , the inter- 
ference of several such terms gives rise to spectral beating 
patterns that underlie shell effects in nuclei, quantum 
dots, or metallic clusters 0. Indeed, as follows from the 
analysis performed by Balian and Bloch [|| , the inclusion 
of just two simplest periodic orbits in a spheroidal cav- 
ity explains the essentials of the shell structure in these 
systems. 

Majority of semiclassical studies on the level-density 
fluctuations was performed for hard-wall systems — two- 
dimensional billiards or three-dimensional cavities 0, Q . 
In these systems, the calculation is considerably simpli- 
fied since each individual orbit exists with easily pre- 
dictable properties for all energies of the particle bounc- 
ing between the walls and contributes by a well-defined 
term to the single-particle level density. Nevertheless, 
the influence of periodic orbits is equally important also 
in systems with "soft" potentials, where the orbit anal- 
ysis is much more involved. This is also the case of 
IBM, where the classical limit for I = describes two- 
dimensional motions within a bounded (for each finite 
E) range of quadrupole deformation parameters, gov- 
erned by a Hamiltonian containing specific kinetic and 
potential terms [1 fiol [Til IT^ . 

The purpose of the present part of our contribution is 
to show that the shell effects and the IBM level bunching 
phenomenon are indeed of similar nature, both originat- 
ing in some particular features of classical periodic tra- 
jectories. Nevertheless, our reasoning does not only point 
to ordinary beating patterns, known from hard-wall sys- 
tems, but makes use of two concepts that in the context 



2 



of nuclear models are somewhat less usual. The first one 
relies on bifurcations of periodic orbits > the second on 
monodromy in classical and quantum integrable systems 
[l4|. Both these effects lead to singular contributions 
to semiclassical trace formulas, that provide a simplified 
description of the level-density fluctuations. Therefore, 
they can be potentially linked to anomalous effects in 
quantal spectra, such as the level bunching at E w 0. 

The paper is organized as follows: In Section |n] we re- 
view the construction of the classical limit of the IBM 
Hamiltonian under study and describe basic features of 
the resulting classical dynamics. Subsection IIII Al briefly 
recapitulates the Berry- Tabor trace formula and the role 
of singular orbits and bifurcations in the semiclassical 
theory of quantal spectra. Numerical analysis of orbits 
with I — 0, presented in Sec. 1111 Bl shows that in the 
£«0 region our system passes through a robust struc- 
tural change of classical dynamics. This change is corre- 
lated with the occurrence of a singular bundle of E — 
trajectories and triggers multiple bifurcations of orbits 
in the region E > 0. The relation of these findings to 
the concept of monodromy is discussed in Section IIVI 
Finally, Section [3 contains concluding remarks. 



II. CLASSICAL HAMILTONIAN 

The classical limit of the IBM can be obtained via the 
well-known procedure, elaborated in detail by Hatch and 
Levit §j, and by Alhassid and Whelan [H 111 El- The 
procedure makes use of Glauber coherent states \a) cx 
exp (a: s st + a^djJlO) with complex time-dependent 
coefficients a = {a s , a M } M= _2,....+2 which define a set 
of 12 classical-like variables (both coordinates and mo- 
menta). The equations of motion for a are derived from 
the time-dependent variational principle, which results 
in the Hamilton function given by the coherent-state av- 
erage H c \(r]]a) = (a\H(rj)\a). This function and analo- 
gous counterparts of other operators can be obtained by 
substitutions s,d^ i— ► ct s ,a^ and s',dl i— > a*, a* in the 
respective normal-ordered quantal expressions. 

Since Glauber coherent states do not fix the total 
number of bosons, an additional constraint must be re- 
quired, namely (a|iV|a) = |a s | 2 -I- ^ M \ a v\ 2 = N. This 
(plus an arbitrary choice of the overall phase, a s — 
N — ^2 \ a n\ 2 ) reduces the number of relevant degrees 

of freedom from six to five. Naturally, the classicality of 
coherent states becomes more and more pronounced as 
N increases and the fully classical limit is obtained in the 
N — > oo limit. To prevent divergence of the correspond- 
ing averages, one has to scale all operators according to 
their order [see the factors in Hamiltonian JJJ] and 
to absorb the respective factors into the definition of ct's. 
This leads to the substitution a,, i— ► a,, = —M? while, 

simultaneously, the ™ factors drop out. 

Final expressions for the classical-limit observables are 
obtained after the identification of real coordinates q^ 
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FIG. 1: Potential energy term of Hamiltonian © for three 
values of parameter rj. The lowermost (77 = 0) and uppermost 
(77 = 1) curves correspond to 0(6) and U(5) limits, respec- 
tively. The middle curve (77 = 0.6) represents an intermediate 
case, for which the accessible range of radii (for the three given 
energies) is shown by the dashed lines. "Negative radii" ex- 
press the (/3,7) — > (/3,7 + 180°) = (— transformation and 
are included just to emphasize the rotational symmetry. 



and momenta p M via relations v25/j = (— ) M( Z-^ + i-Pn 
and V2a* = — (— The coordinates q^ are 
associated with the geometric variables describing an in- 
stantenous quadrupole deformation of the nucleus and 
its orientation in the laboratory frame. Due to the fixed 
boson number average, the motion is constrained by the 
condition 



(2) 



to the interior of a sphere in the 10-dimensional phase 
space. 

The calculation of classical observables is substantially 
simplified for zero angular momentum, I = [T^ . Il5j . In 
this case, the intrinsic frame connected with the ellipsoid 
of deformation remains at rest and one can fix Re qo = x, 
and Req +2 = Reg_2 = y/V% (while q±x = Img±2 = 
Imgo = 0). The I = classical limit of Hamiltonian (JTJ 
reads as 



Hi = 7 ]7r 2 + (l-r ? )^V -i 



2 1 5v -f + (i-v)P\ 



(3) 

where 2 — x 2 +y 2 is the squared radius in the qo x \/2q±2 
plane (the polar angle denoted as 7) and ir 2 the squared 
length of the associated vector of momenta: 



7T = 7T X + TT y = ftp 



(4) 



Hamiltonian @ can be thought to describe planar mo- 
tions of a particle with the position-dependent kinetic 
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energy T c \ in potential V c \, which is for 77 = 0, 0.6 and 1 
shown in Figure ^ While for 77 < rj c = |, the potential 
has the "Mexican hat" (or "champagne bottle") form, 
for 77 > it is just a well with minimum at [3 = 0. To 
emphasize the rotational symmetry in the x x y plane, 
we show both positive and "negative" domains of (3 (the 
latter corresponding to the rotation by angle 180°). As 
follows from the form of the potential and from Eq. J2Jl, 
the radius must satisfy 



P G [Anin, U C [0, y/2] 



and the total energy 



E G [-E m in, £max] C [— 1, +1] , 



(5) 



(6) 



En 







(5»7-4) 2 
16(1-7?) 



for 7] < 
for 7] > 



E u 



Remind that throughout this paper the energy is always 
expressed in units of the scaling constant a, see Eq. flip, 
so it is formally dimcnsionless. 

Note that polar coordinates [5 and 7 of a; and y can be 
immediately associated with Bohr geometric variables, 
but in this case the deformation parameter would be re- 
stricted to {3 G [0, v2). To obtain (3 S [0, 00), as is usual 
in nuclear structure, the coordinate plane must be radi- 
ally stretched IJj according to /3 1— ► (3 = —j£==. Ex- 
pressed in new polar coordinates (f3, 7) and the associ- 
ated momenta (7^73, 7r 7 ), Hamiltonian JSJ transforms into 
the following form: 



1 r 
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r) + (4 - 3t7)/3 s 



(l + /3 2 ) 2 ^+^ 

(5r?-4)/3 2 +?7/3 4 
(1 + /3 2 ) 2 



(7) 



Here, the upper physical limit of energy, E max , is reached 
when the motion becomes infinite. In the following, we 
will use the classical limit in the form of Eq. ©, i.e., with 
f3 restricted to the finite interval (JSJ. 

It is immediately apparent that the Hamiltonian in 
Eqs. (0 and Q) is "7-soft", invariant under rotations 
about the origin, so it conserves the "angular momen- 
tum" 



XTTy - yir x . 



(8) 



Thus, since the number of degrees of freedom / = 2, 
the system must be integrable. This is in agreement 
with the arguments explaining the integrability of the 
[0(6)— U(5)]dO(5) Hamiltonians with arbitrary angu- 
lar momenta, as outlined in Part I 1], since Eq. © 
is closely related to the integral of motion (72[0(5)] = 
i(i • L) +2(f 3 • f 3 ). Indeed, for I = the classical limit 
of the 0(5) Casimir invariant reads as 0: 



C 2 [0(5)] c i 



(9) 
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FIG. 2: Poincare phase-space sections for Hamiltonian Q 
with 77 = 0.6 at the three given values of energy (panels a-c). 
The sections show a finite number of crossings of 10 (a), 30 
(b), and 50 (c) trajectories with the plane /3 x 7173 for arbi- 
trary 7. 



Note that since "angular momentum" ljUJi does not corre- 
spond to ordinary 0(2) algebra of two-dimensional rota- 
tions, its quantization yields eigenvalues v(v + 3), where 
for I — the seniority takes values v — 0, 3, 6, . . ., in 
contrast to the m 2 formula with m = 0,±1,±2, . . . cor- 
responding to 0(2). Nevertheless, we realize that each 
value of C?2[0(5)] c i is associated with both signs of 7r 7 , 
i.e., with two opposite orientations of the motion in 7- 
direction. This intrinsic "degeneracy" (which does not 
affect physical results in the quantum case) will become 
important in Sec. IIVI 

The integrability of Hamiltonian © is illustrated in 
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Figure where we show Poincare phase-space sections 
for 77 = 0.6 at three different energies, (a) E = —0.1, 
(b) E — 0, and (c) E = 0.1. Each of the panels rep- 
resent passages of 10-50 randomly selected trajectories 
with the given energy E through the (3 x irp plane in the 
4-dimensional phase space. Due to the rotational sym- 
metry, the plane can have an arbitrary orientation in the 
x x y frame and the pattern of sections must be symmet- 
ric under the reflection of the j3 axis (we show both (3 > 
and "/? < 0" halves). 

All sections in Fig.|21demonstrate fully regular dynam- 
ics, in agreement with the integrability of our system. 
As can be anticipated from Fig. ^ the E < motions 
in panel (a) must be confined inside the annular region 
[3 S [/3 m in, /Smax], while the E > trajectories in panel (c) 
already range over the full disc (3 € [0,/3 max ] (the values 
/3 m in and /3 m ax depend on energy). Panel (b) shows just 
the singular E = situation, when the central inacces- 
sible (for E < 0) disc shrinks into a single point (which 
can be reached in infinite time). We will return to this 
case in Sec. I1VI 

Let us stress that the Poincare sections in Fig. |21 sepa- 
rate trajectories with different energies, but mix together 
those with various values of the angular momentum 7r 7 . 
Indeed, the outermost curves in all panels represent pure 
/3-vibrations with 7r 7 = 0, while the central points (not 
shown) correspond to "spinning" only in the 7-direction 
with up = 0. The other trajectories correspond to various 
mixtures of (3- and 7- vibrations. An interesting attribute 
of these intermediate cases is the spread A/3 of each indi- 
vidual trajectory in the /3-direction, which can be deter- 
mined as the difference between radii corresponding to 
the outermost and innermost points. This will be used 
in Sec. IHII to classify quasiperiodic orbits in our system. 



III. PERIODIC ORBITS 



A. Berry- Tabor formula, singular orbits, and 
bifurcations 



Semiclassical analyses of quantal spectra are performed 
in the framework of so-called trace formulas which rep- 
resent the fluctuating part ga(E) of the level density 
£>tot = Qsm + Qa in terms of purely classical quantities as- 
sociated with periodic orbits, while the complementary 
smooth part g sm (E) is determined just as the classical 
phase-space volume available at a given energy E in units 
of 7i/ (where / is the number of system's degrees of free- 
dom). The best known expression, derived by Gutzwiller 
P|, was obtained under the assumption that individual 
periodic orbits are isolated, which is not satisfied for in- 
tegrable systems. In this case, periodic orbits come in 
continuous families characterized by arbitrary shifts of 
initial angles if the motion is described in the action- 
angle variables |l3j . An adequate semiclassical approach 
to the level density of integrable systems was developed 
by Berry and Tabor . 



In the following, we consider a two-dimensional inte- 
grable system, / = 2, which applies in our case of Hamil- 
tonian Q • In the action-angle representation the Hamil- 
tonian depends only on actions, H = H(Ii, I2), and the 
angles evolve according to Q\ = uii and 82 — u>2- All 
trajectories represent folded rotations on various tori de- 
termined by I = (7i, J2). Let us note that in many inte- 
grable systems, including ours, the action-angle variables 
can only be introduced locally 0| (cf. Sec. llVj l. There- 
fore, it is not possible to write down analytic expressions 
for the corresponding canonical transformation from nor- 
mal coordinates and momenta. 

Any primitive periodic orbit on a given torus can be 
characterized by a pair of coprime integers (jii,^) = M 
such that the ratio of angular frequencies 1Z = ^ co- 
incides with the rational number j^. The Berry- Tabor 
formula for the fluctuating part of the quantal state den- 
sity [3 then reads as 



g&(E) 



1 00 

~Trh 



Tit 



rSjt(E) 



(10) 



where the sum runs over all repetitions r of all primitive 
orbits a with period Tn — 27rpi = 
jjj, and action 



Maslov index vp 



S fl (E) = 2irl ■ (1= [71(3/3 + ^ 7] dt 



T- 



(11) 



The meaning of the function g'^ in Eq. (|10|) will be ex- 
plained later. 

Expression l|ll|) . which in the general case integrates 
the scalar product of momentum and velocity over the 
specific periodic trajectory fl, has a particularly sim- 
ple form for billiards (or cavities), where one can write 
Sp = 2ETp = pLp with p — mv denoting the ordinary 
momentum and Lp the length of the given orbit. For 
"soft" systems, the dependence of Sp on energy is nonlin- 
ear and frequencies of individual cosine terms in Eq. (|10|) 
vary with E. Since in the latter case each oscillatory 
term in the Berry- Tabor formula contains also a nontriv- 
ial energy dependence of the amplitude, the semiclassical 
analysis of spectra in such cases is certainly much less in- 
tuitive than in the hard- wall systems. 

In general, there may exist singular orbits with diverg- 
ing contributions to the Berry- Tabor formula. This hap- 
pens if either the period of the given orbit grows to in- 
finity, Tp — > 00, or if the denominator of the prefactor 
in Eq. (|10fl vanishes, g'^ — » 0. The former case applies 
to the motions that for some energy become infinitely 
slow at a certain point, which can be associated with 
an unstable equilibrium of the system. We know from 
the discussion in Sec. ITll that our system contains such a 
point, namely the central maximum of the potential in 
Eq. 01 at (3 = for 77 < |. For trajectories with E = 0, 
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this maximum can only be reached in asymptotic times 
because the force vanishes there. Among the trajectories 
passing this point there are also various periodic orbits, 
whose contribution to Eq. (|ll)fl must diverge at E = be- 
cause of the period tending to infinity. This is essentially 
the classical mechanism responsible for the bunching of 
quantum levels in the region E as 0, see Fig. 1 in Part I 
It shows that the bunching pattern is not just a 
finite-iV quantum fluctuation, but a robust effect deeply 
ingrained in the classical limit of the system. Theoretical 
foundations underlying the existence of the singular class 
of trajectories and another approach to understand their 
influence on the quantum spectrum will be discussed in 
Sec.lTvl 

The second possible source of infinite contributions to 
the Berry- Tabor formula I|1U|) is connected with the cases 

when g' E (h) = ^jf-(h) = 0. The function g E {h) is 

determined from the implicit equation H(Ii,I 2 — 
9e) — E>, which after differentiation and the use of Hamil- 
ton equations yields 



' dh 



, 



(12) 



so that g' E = — ^ — —7Z. In other words, the function 
gE matches possible pairs of actions (ii, J2), i.e., selects 
the tori / relevant at a given energy, and its first deriva- 
tive determines the corresponding frequency ratios. If 
— g' E is rational, for a selected torus, the associated or- 
bit is periodic and contributes to Eq. I|1U|) . The second 
derivative g' E measures the change of 1Z as one steps to 
the tori in an infinitesimal vicinity of /. If g E ^ 0, the 
periodic orbit /2 on the torus / does not survive the tran- 
sition to / + SI. If, however, g E = 0, a family of periodic 
orbits with the same frequency ratio 1Z exists in neigh- 
boring tori, which results in diverging contribution to the 
Berry- Tabor formula (|10fl . 

Note that the Gutzwiller formula [Bj], which is valid 
in nonintegrable systems with isolated periodic orbits, 
is formally similar to Eq. I|10[l . but with the prefactor 
denominator replaced by ^/det[(A/ p ) r — 1], where M p 
stands for the so-called monodromy matrix of a given 
primitive periodic orbit p 0, Q . This matrix describes 
the stability of orbit p in terms of linearized deviations 
from the given phase-space trajectory under a perpendic- 
ular perturbation of its initial point. Thus (M p ) r — 1 rep- 
resents the deviation from the perturbed phase-space po- 
sition after r repetitions. If one (or more) of the eigenval- 
ues of this matrix is equal to zero, i.e., if det[(M p ) r — 1] = 
0, there exists at least one direction in the phase space 
in which any deviation from the given orbit r ■ p results 
in another periodic orbit. The new orbits are detached 
from the primitive orbit p as its period r-tupling clones. 
Consequently, p is not isolated and the corresponding- 
term in the Gutzwiller formula diverges. This situation 
is analogous to the one with g E — 0, as described above. 

Both the above singular cases correspond to the same 
general phenomenon, called bifurcation |18| . In Hamil- 



ton systems of classical mechanics, bifurcations repre- 
sent branching of periodic orbits at some critical values 
of energy or other parameters 0]. While periodic or- 
bits existing below and above the given bifurcation en- 
ergy £"b are isolated, either in the sense of g E ^ or 
det[(Af p ) r — 1] 7^ 0, at E = E\> two or more orbits merge 
together in the way described above, giving rise to zero 
denominators of the respective semiclassical level-density 
formulas. At the bifurcation energies E^, the Berry- 
Tabor or Gutzwiller formulas do not represent correct 
approximations of the fluctuating level density. Improved 
semiclassical methods were developed to treat these sit- 
uations |l9j| . Intuitively one expects an enhancement of 
the level-density oscillations at the bifurcation points. In 
the following subsection we will show that in our system 
numerous bifurcations of periodic orbits take place in the 
energy range E > 0. 



B. Numerical results 

We have performed a numerical analysis of classical 
motions corresponding to Hamiltonian in the inter- 
val of energies E e [—0.1, +0.3] using a sample of about 
50000 generated orbits. Individual trajectories were cal- 
culated with initial positions and momenta chosen ran- 
domly within the phase-space region accessible at a given 
energy and classified by the ratio 



R=^ 

Tf3 {uj 



{up) 



(13) 



7/ 



of periods T 7 and Tp associated with oscillations in both 
7 and /3 directions, respectively. Since w 7 = 7 and the 
angular velocity u)p connected with /3-vibrations are both 
time dependent, one has to use the corresponding average 
angular frequencies per period, (w 7 ) = -f- and (ujp) = 

. Their inverse ratio coincides with R and is analogous 



to the above-discussed ratio 1Z 



of frequencies in the 



action-angle variables. In particular, rational values R = 
j& correspond to periodic orbits with period T^j^ = 
UpTp = n 1 T 1 . 

Examples of periodic orbits with various rational val- 
ues of the ratio l|13f) are shown in Figure |3| The rational 
fraction R — classifying the given orbit has visual 
meaning as the number of outer return points of the j3- 
vibration over the number of rotations in the 7-direction 
needed to close the orbit. Thus, for instance, the 5/2 
orbits look like stars with 5 outer "points" that close in 2 
rotations, while the 5/1 stars are similar, but close only 
in 1 rotation. The outer and in some panels also the in- 
ner circles in Fig.|21demarcate the energetically accessible 
areas (3 £ [/3 m in,/3max] in the x xy plane. As discussed in 
Sec. [n] this area is a disc for E > 0, an annular ring for 
E < 0, and a disc minus the central point for E = 0. We 
see that although the orbits in Fig.[3]do not just trivially 
bounce between the outer (and inner) limits, as in the 
case of circular or annular infinite sharp wells, they still 
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FIG. 3: Examples of various periodic orbits at different ener- 
gies for Hamiltonian © with r\ — 0.6 and their classification 
by rational fractions R = — . 

resemble to a large extent the trajectories in these simple 
systems pflj . 

Periodic orbits form a dense subset of all allowed mo- 
tions and we therefore need a more complete picture. 
A histogram showing the occurrence of trajectories with 
arbitrary (rational or irrational) values of the frequency 
ratio R within the whole sample of trajectories with 
E G [—0.1, +0.3] is presented in Figure^] For each value 
of energy within the given range (the energy step was 
chosen AE = 0.01), the sample contains iVtot = 1200 tra- 
jectories and Fig. 0] depicts their distribution (numbers 
Ntr of trajectories) into bins of size AR = 0.01 along the 
i?-axis. 

The structure shown in Fig. 0] discloses rather inter- 
esting features of classical motions. For each energy, the 
orbits occur within a band 

Re[R min (E),R m , x (E)] (14) 

of allowed frequency ratios. The lower bound R m - ln (E) 
gradually decreases with increasing energy for E < 0, 
but it is constant, i? m i n = 2, for E > 0. Because of 
the limited precision inherent in our generated sample of 
trajectories (with nonzero values of the bin size AR and 
energy step AE) we cannot resolve whether the i? m j n (i?) 
dependence is discontinuous or just nonanalytic at E = 
0. In any case, the curve reaches the minimal value at 
this point. On the other hand, for energy approaching 
the minimum Eq of the potential in Eq. Q (this energy 
is below the range displayed in Fig. we must have 
^min - > co. The upper bound of interval (|14[) is also 
a decreasing function of energy which passes the value 
fl max = 4 at E = 0. The decrease of R max (E) for E < 
(and partly also just above E = 0) is so steep that it 
cannot be resolved with the present energy step, but we 
assume that it is a smooth curve. It is obvious that a very 
narrow energy interval around the point E ~ carries the 
most substantial changes in the spectrum of orbits, where 
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R 

FIG. 4: Frequency of occurrence of trajectories with different 
ratios R for Hamiltonian ©, r\ = 0.6, in the reference sample 
of trajectories with different energies. For each value of E 
(step AE = 0.01) there was 1200 generated trajectories and 
the histogram (see the inset) shows their distribution in R 
(the bin width AR = 0.01). The main diagram (contour plot 
of the logarithmic histogram) depicts the band structure of 
allowed R values, see Eq. 11411 and below. 



the trajectories pass between both negative and positive 
energy regions just through a bottleneck of values R G 
[3, 4). While for E < the orbits look similar to those in 
the 0(6) limit, for E > they already resemble the U(5) 
limit. 

The behavior demonstrated in Fig. 0] can be qualita- 
tively understood from the change of the energetically 
accessible x x y area around E w 0. The form of an- 
nular ring, valid for E < 0, does not support trajecto- 
ries with R < 3 since these have to traverse through the 
central region. Consequently, these trajectories can only 
exist for E > —e, where e ~ 0.03. On the other hand, 
the central reflecting disc is needed for trajectories with 
R > 4, which therefore appear only for E < 0. One can 
say that the R G [4,oo) trajectories, which are "bounc- 
ing" between inner and outer circles inside the annular 
region for E < 0, transform to the straight R G [2,3] 
trajectories at E « where the central disc gradually 
disappears and the accessible domain of deformation pa- 
rameters becomes simply connected. Note that the ra- 
pidity of changes of classical motions around zero energy 
is connected with the fact that for E — > 0_ the radius of 
the central disc converges to zero with a rate increasing 
to infinity (/3 m i n oc \/—E), as directly follows from the 
form of the potential V^i close to the (3 = maximum. 

Figure shows the relative frequency of occurrence of 
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(b) 



0.1 0.2 0.3 -0.1 

energy 

FIG. 5: Relative frequency of occurrence of several types of 
quasiperiodic orbits (see Fig. [3J for r\ — 0.6 as a function of 
energy. The scale on the vertical axis depends on the value 
A7 in Eq. 115L here equal to 5°. 



several types of periodic orbits from Fig. [3] in our gen- 
erated sample as a function of energy. The curves in 
Fig. [5] can be basically understood as energy cuts of the 
function N tr in Fig. 0] at the respective rational values 
R = — of the fraction (II 31) . but with a variable preci- 
sion AR. More specifically, Fig. 03 presents the relative 
fraction of all generated (at each energy) trajectories sat- 
isfying the condition that the ^t/jth outer reflection after 
fi-y revelations is shifted from the 1st outer reflection by 
an angle not exceeding (in absolute value) the selected 
precision A7 = 5°. This leads to the condition 



Hi 



R 



< 



A7 

27T/i 7 



(15) 



Let us stress that the use of a smaller value of A7 de- 
creases the yield of trajectories — implying a prolongation 
of the computation time — but does not change (as we 
checked for A7 = 1°) the shape of dependences in Fig. E] 

The most common type of behavior shown in Fig. 
indicates that for many orbits the relative frequency of 
occurrence sharply culminates at a certain energy, just 
before this orbit totally disappears from the system. The 
sharpest peak of this kind is observed at E = 0.02 for 
the 3/1 orbits, but there are also other well pronounced 
peaks, like the 5/2, 11/4, or 13/5 ones, and many others. 
All these maxima appear at positive energies and one 
can trace their origin to the ridge of values N tr visible 
for E > at the upper bound R max (E) in Fig. H (see the 
inset). The peak at E = (the 4/1 "crosses") and also 
the one at E = 0.02 (the 3/1 "Mercedes-Benz stars", see 
Fig. |3J) are located just on the upper edge of the major 
£»0 level bunching pattern in Fig. 1 of Part I 0. 

Special attention should be paid to the 4/1 orbits that 
in our system take two different forms: For E < they 
exist as stars, shown in the third uppermost panel of 
Fig. |3 but at E w they can also look like crosses, see 
the second panel. (In fact, the latter case exemplifies the 
above-discussed critical E = periodic trajectories with 
infinite period, as will be further elaborated in Sec. II VP 
The contributions of these forms to the dependence in 
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0.3 
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-0.1 



3/1 i — i — 

8/3 — x~ 

12/5 >--*- 

7/3 i a 




0.6 



1.2 




radial width 



FIG. 6: The radial width A/3 of individual periodic orbits — 
as a function of energy (»y = 0.6). Horizontal bars at each 
energy demarcate intervals of A/3 where trajectories are de- 
tected within our sample, while points connected by curves 
represent statistical averages within each interval. The bifur- 
cations of orbits in panel (a) explain the respective peaks in 
Fig.El 



Fig. Ela) can be decomposed into a constant step-like 
function equal to zero for E > ("stars") and a sharp 
peak at E = ( "crosses" ) . 

We also see in Fig.|5Ja) that the 2/1 orbits, which pass 
via the central maximum of the potential and correspond 
to the E > edge i? m i n in Fig.0J exhibit a different type 
of energy dependence than the others. The frequency of 
occurrence of these orbits is zero at E < and gradually 
increases (if neglecting fluctuations) with energy E > 0. 

It is not difficult to show that the peaks in Fig. [5] cor- 
responding to the R £ (2,3] orbits are connected with 
bifurcations. To this end, we first characterize individual 
orbits by the radial width, defined as the difference A/3 
between the outer and inner radii (see the end of Sec. lTT|l . 
The values of A/3 associated with various orbit types are 
shown in Figure where the horizontal bars demarcate 
intervals of the A/3 values which are populated (for en- 
ergy given on the vertical axis) by some trajectories in 
our sample. The curves (used just to lead the eye) con- 
nect points that represent arithmetic averages of A/3 in 
neighboring intervals. 

The three panels of Fig. [fj] collect three types of qual- 
itatively different behaviors: (a) For R £ (2,3], the pop- 
ulated domain of A/3 consists of two separate branches 
(see, e.g., the two 7/3 orbits in Fig. [3J that merge at a 
certain energy E-b(R) > 0, which can be determined from 
the condition i? ma x(-£'b) = R- Similarly, the lower end- 
point energy (the termination of the longer branch of the 
respective curve) follows from the R m i n (E) bound, (b) 
For R £ (3, 4], the domains consist of only one band that 
shifts to larger A/3 values as the energy increases and 
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terminates slightly above E — 0. The upper endpoint 
energies can again be determined from R max {E), but as 
discussed above, this dependence is so steep in the given 
range of R that all endpoint energies in Fig. HJb) fall 
into the narrow interval E £ [0,0.01]. Lower endpoint 
energies again follow from R m - m (E). A special case of 
this kind is the 4/1 orbit with the two above-discussed 
incarnations (see Fig.JSJ: the respective A/3 value at the 
endpoint E — is apparently deviated from the direction 
followed for E < 0. (c) For R > 4, the A/3 bands look 
similar as in case (b), but terminate at energies just be- 
low E = 0, supposedly following the steep E < branch 
of the curve R m ax(E). 

It becomes apparent that the R € (2, 3] peaks in Fig.0 
arise due to the merge of two different branches of A/3 
values, as shown in Fig.EJa). An increased frequency of 
occurrence of the orbit just before the endpoint is related 
to the flatness of the respective curve close to its maxi- 
mum (a larger number of trajectories is concentrated in 
a smaller energy interval). At the endpoint (A(3 b ,E b ) of 
each of the curves in Fig. |5f a) the respective type of pe- 
riodic orbit bifurcates, having infinitely close neighbors 
with different radial widths, and thus yields = 0, as 
discussed in Sec. IIII Al The Berr y- T abor formula Hl()(l 
cannot be applied at these points [l£|- In our case, the 
bifurcations seem to be of the pitchfork type, when two 
stable orbits join and produce an unstable one fl8| . Un- 
fortunately, the unstable orbits are not accessible to nu- 
merical studies, so they are not seen in Fig. This 
problem may be further investigated analytically. 

It follows from the above discussion that the bifur- 
cations are connected only with the region of positive 
energies (there is a ridge of N tT values, apparent in the 
inset of Fig. 0] which is located solely at the R max (E) 
edge with E > 0) . This implies that divergences of the 
Berry- Tabor formula associated with bifurcations are not 
directly relevant in the explanation of the main level- 
bunching pattern in Fig. 1 of Part I (except perhaps 
the 3/1 case with E b = 0.02). Bifurcation energies for 
low-period orbits are not even correlated with the sec- 
ondary, less pronounced bunchings of levels, observed in 
the region E > [J]. Therefore, it seems that the pres- 
ence of various orbits in the same energy range and an 
interplay of their bifurcations result in interferences that 
wash out contributions of individual orbits. 

On the other hand, highly organized behavior of levels 
at E perfectly coincides with the predicted existence 
of a singular torus of orbits with infinite period at zero 
energy, and also with the observed abrupt redistribution 
of the spectrum of orbits in a narrow vicinity of this en- 
ergy. 



IV. CLASSICAL AND QUANTUM 
MONODROMY 

The anomalous E = bundle of orbits with infinite pe- 
riod, discussed in Sec. IIII Al is related to a more general 



phenomenon, called monodromy. Classical monodromy 
in integrable Hamilton systems can be briefly introduced 
as the impossibility to define global action-angle variables 
due to the existence of a singular, so called "pinched" 
torus 01- The name MovoSpofitt ("once around") orig- 
inates from a property similar to that of the Mobius 
strip: if one follows a closed loop in the space of regular 
tori around the singular torus and — loosely speaking — 
redefines the coordinate system on the consecutive tori 
continuously on the way along the loop, one returns back 
to the starting torus with a coordinate system that differs 
from the initial one. 

Classical monodromy affects the quantum counterpart 
of the system via the Einstein-Brillouin-Kramers (EBK) 
quantization rules 0. Of ■ It turns out that quantum mon- 
odromy can be seen as a point defect in the lattice of 
quantum numbers corresponding to a complete set of 
commuting operators. This defect results in a transfor- 
mation of the elementary quantum cell when a closed 
loop is completed around the singular point, in analogy 
with the above feature of phase-space tori. An overview 
of the mathematical background and various examples of 
monodromy can be found in Ref. |2l| . 

Soon after its discovery in 1980 it became clear 
that monodromy substantially affects global features of 
numerous integrable systems, which might previously be 
considered as too trivial for detailed analyses. The sim- 
plest system that exhibits monodromy is the spherical 
pendulum — particle moving on a sphere in a gravitational 
field. It can be shown 123 that the phase -space torus 
passing the unstable equilibrium position at the north 
pole, with the particle energy exactly equal to the criti- 
cal value E m needed to reach that point, is pinched, i.e., 
one of its basic circles is contracted to a single point (with 
appropriate initial conditions the particle is at rest). As 
a consequence, the lattice of quantum states, character- 
ized by quantum numbers enumerating energy E and the 
projection L z of angular momentum, has a point defect 
at (E, L z ) = (E m , 0). It was found that closely related to 
this simple observation is the realization of monodromy 
in vibrational and rotational spectra of some molecules 

mm. 

Other examples of monodromy can be found in the 
following s yste ms: particle in quartic, sextic, and decatic 
potentials |24| , hydrogen atom in orthogonal electric and 
magnetic fields 1251 . systems of two or three coupled an- 
gular momenta .26], particle bouncing between walls in 
a prolate elliptic cavity |27| or moving in a two-center 
attractive potential [28j. Like in the spherical pendulum, 
monodromy in several of the latter systems is connected 
with the trajectories passing with the critical energy via 
the point of an unstable equilibrium [2ll I24I EsT ] . We al- 
ready know that a similar point, namely the top of the 
central maximum at (3 = 0, can be found in our classi- 
cal Hamiltonian (J3J for rj < r\ c = |, with E m — being 
the minimal energy needed to pass through this point. 
In fact, our Hamiltonian for rj < rj c is identical with the 
"champagne-bottle" Hamiltonian of Ref. [24| , except the 
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20 40 20 40 20 40 



Vv(v + 3) 

FIG. 7: The lattice of Z = eigenstates of Hamiltonian Q 
with three given values of r\ (panels a-c) and N = 40 in the 
plane E x v/ 'v(v + 3) (where the seniority v = 0,3,6,...). 
Lines connect states with the same radial quantum numbers 
rip. The singular torus (E, v) — (0, 0) is located in the center 
of the semicircle in panel (b). 

position-dependent kinetic term in Eq. @, which how- 
ever does not affect the presence of monodromy with the 
central point (E,ir^) = (0,0). 

Figured shows the lattice of I = eigenstates of quan- 
tum Hamiltonian with the number of bosons N = 40 
in the plane where the vertical axis represents energy E 
and the horizontal axis the momentum 

tt 7 = y/v(y + 3) = 3y/v(v + 1) . (16) 

Since for zero angular momentum the seniority takes val- 
ues equal to multiples of 3, we defined above also the 
"reduced" seniority quantum number v = ^ = 0, 1, 2, . . . 
The three panels in Fig. correspond to various values 
of the control parameter: (a) rj = 0, the 0(6) case, (b) 
i] = 0.6, a transitional case, and (c) 77 = 1, the U(5) 
case. Individual states (marked by dots) can be directly 
related to level energies at the respective values of 77 in 
Fig. 1 of Part I p] , which collects all states with different 
seniorities for the same boson number as here. 

Sorting of states according to seniority in Fig. helps 
to identify the values of v that are involved in level bunch- 
ings at different points 77. For instance, one immediately 
sees that the clustering of levels across the whole spec- 
trum in the U(5) limit (panel c) is due to the multiple 
degeneracy of states with even or odd values of v, that 
correspond to the same value of the U(5) quantum num- 
ber rid (even or odd, respectively). For the highest states, 
this degeneracy remains approximately valid across the 



whole interval 77 S [0, 1], see panels (a)-(c). The seniority 
deconvolution of the spectrum for 77 < | is exemplified 
by the 77 = 0.6 case in panel (b). We observe here that 
levels with all values of v become nearly degenerate in 
the region around zero energy, which is a clear signature 
of the E ss bunching pattern Q . 

The lattices in Fig. represent quantum energy- 
momentum maps [2l| of the classical phase space, with 
each dot being an image of a classical torus of trajectories 
that survived the semiclassical EBK quantization [l7j . 
This is given by E = 2wh(rii + ^j) , where Ii (with i = 1,2) 
are quantized actions and V{ the respective Maslov indices 
£j . The EBK tori should be determined by two quantum 
numbers ri\ and 77,2 > whose integer values increase by one. 
Good candidates for these numbers are the reduced se- 
niority v (connecting vertical columns of points in Fig. 
and the radial quantum number rip — 0,1,2,..., that 
enumerates states with a fixed v according to energy (in 
Fig. the constant-n^ states are connected by lines). 
For 77 = 0, the radial quantum number is related to cr, 
which corresponds to the 0(6) Casimir invariant [2|, and 
the pair (np,v) represents the appropriate choice of the 
EBK quantum numbers. In the U(5)-like shown 
below, yet an alternative pair of quantum numbers needs 
to be defined. 

It follows from Eq. that the I = classical limit of 
the 77 = 1 Hamiltonian (0) is identical with an isotropic 
two-dimensional harmonic oscillator. Indeed, for the sub- 
set of states with the U(5) quantum number rid equal to 
multiples of 3 (in this case rid = 2np + v) the U(5) lat- 
tice coincides with the 2D-oscillator lattice of states (the 
energy in the latter case being enumerated by the oscil- 
lator quantum number n a = 2n r + m, where n r and m 
stand for ordinary radial and angular-momentum quan- 
tum numbers, respectively). In the entire U(5) lattice, 
however, majority of states is located in "interstitial" 
positions with rid 7^ 3/c; this is because the underly- 
ing "angular-momentum" algebra differs from the ordi- 
nary 0(2). Apart from rip and v, all U(5) states can 
be labeled by a pair of oscillator-like quantum numbers 
Til = + v and 77,2 = np + 2tj. States with constant val- 
ues ni — 0,1,2,... form upwards inclined rows of dots 
in Fig. 0c), while the ri2 = 0, 1, 2, . . . quantum number 
connects states in the downwards inclined rows. The ri\ 
chains are clearly apparent also in both remaining panels 
(a) and (b) of Fig. 

For the purpose of the semiclassical analysis, the lat- 
tices in Fig. must be extended to cover both positive 
and negative 7r 7 . Remind from Sec. ^ that although 
the physical quantum states can be represented by non- 
negative values of 7r 7 , the intrinsic degeneracy of classical 
motions in both 7-directions results in mirror imaging 
of all states with v > into the 7r 7 < half-plane (to 
guarantee a smooth continuation of quantum numbers, 
we assign values v = — 1, — 2, ...to these "twin" states). 
In absence of monodromy, one must be able to engage 
all states in the extended lattice into a "crystal" grid of 
continuous and smooth lines, corresponding to constant 
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values of two compatible global quantum numbers, with 
"elementary cells" of the grid being topologically equiv- 
alent to squares. From Fig. 0[a) we see that a smooth 
grid, symmetric under the 7r 7 <-> — 7r 7 reflection, can be 
constructed in the 0(6) case, using the pair of generat- 
ing quantum numbers (n/3, v). In the U(5) case (panel c), 
this choice of quantum numbers produces a grid of lines 
that are broken at 7r 7 = 0, but a smooth global grid (a 
diagonal "chessboard") is generated by the pair (711,712). 
The latter structure can be extended to the whole inter- 
val 77 € [|, 1] where the U(5)-like spectrum exists. 

In contrast, quantum monodromy implies the absence 
of a smooth global grid. This is the case of Fig. 0[b), 
where a smooth grid for E < would be generated by 
the pair of quantum numbers (np,v), but for E > by 
the pair (ni, 712). Any attempt to define two global quan- 
tum numbers that behave smoothly in the entire lattice 
for 77 6 (0,1) fails at the point (E,v) = (0,0), which 
represents the singular torus of trajectories and, simul- 
taneously, a "defect" in the quantum lattice of states 
[2HI24I l2q . It is clear that in the transition to the 0(6) 
limit the defect is gradually pushed up to the upper edge 
of the lattice. For the whole interval 77 € (0, 4) the singu- 
lar point indicates the place where the energy-momentum 
map passes between the 0(6) and U(5) types of elemen- 
tary cells — tetragons with ordered {np 1 v) 1 {np,v + 1), 
(np+1, (n^+l, v) vertices, and analogous tetragons 

in (711,71.2), respectively. Note that elementary cells of 
either type cannot be uniquely defined along a closed 
loop around the singular point since after one turn the 
cell gets distorted. This can be illustrated by a graph- 
ical construction in Fig. Ctb) and its mirror image, but 
the rigorous proof would require an infinite density of 
the lattice in the N — ► 00 limit. The last observation 
represents a common quantum signature of monodromy 

s mm m mm. 

It should be stressed that monodromy in the present 
case is not a property of just a single Hamiltonian, but 
characterizes the whole r\ € (0,?y c ) family JIJ of transi- 
tional systems. Since — as shown in Part I [lj — the most 
substantial changes in quantum spectra of these systems 
take place in the E sa region, monodromy seems to 
play the key role in the process of redistribution of indi- 
vidual levels between the 0(6) and U(5) multiplets. Re- 
lated examples exist also in other parametric families of 
Hamiltonians, for instance, in transitions between uncou- 
pled and coupled regimes of two quantum rotators j2|| 
and between Zeeman and Stark limits of the hydrogen 
atom in crossed electric and magnetic fields |2j| . Also in 
these examples, the crossover between the limiting spec- 
tral structures takes place at the point (or in the interval) 
of control parameters and energy where monodromy ex- 
ists. These findings deserve further investigation. 
V. CONCLUSIONS 

In the present part of our work, devoted to the 
[0(6)— U(5)]dO(5) transition of the interacting boson 



model, we have studied the classical limit of Hamilto- 
nian Q with zero angular momentum. Results of the 
analysis of level dynamics, presented in Part I , were 
qualitatively discussed with the aid of the Berry- Tabor 
semiclassical trace formula, which describes fluctuations 
of quantum spectra in integrable systems in terms of fam- 
ilies of periodic orbits existing at various energies. The 
transitional regime was exemplified by the choice of a 
single value of the control parameter, r\ = 0.6. 

Both possible sources of diverging contributions to the 
trace formula l|l()|l . namely, the existence of bifurcating 
(g'g = 0) and singular (Tp — 00) orbits, were identi- 
fied in our system. While bifurcations of periodic orbits 
with ratios between 7- and (3- vibration periods R € (2, 3] 
were shown to exist in the region E > 0, singular orbits 
with 7r 7 = appear at E = 0. The latter finding led 
to the identification of classical monodromy, which on 
the quantum level exhibits itself as a defect in the lattice 
of quantum states located at zero values of energy and 
seniority. This results in the bunching of levels in the 
E s» region Q and underlies the process of redistribu- 
tion of states between 0(6) and U(5) spectral structures, 
i.e., between the (np,v) and (711,712) types of elementary 
cells, and the respective multiplets of levels. 

Also the numerical analysis of periodic and nonperiodic 
classical vibrations disclosed that the most substantial 
changes in the spectrum of allowed ratios R take place 
in a very narrow energy interval around £«0. This in- 
terval represents a kind of demarcation line between the 
0(6) and U(5) types of classical motions. At E w —0.03, 
the inaccessible central disc in the plane of deformation 
parameters becomes sufficiently small to allow for vibra- 
tions with R e (2, 3], and at E = 0, when the disc van- 
ishes, the oscillator-like orbits with R = 2 arise. With 
the energy further growing to positive values, individual 
vibrations with R > 2 eventually disappear in bifurca- 
tions ("annihilations" of two separate A/3 branches of a 
given orbit). 

We believe that results of this analysis will have a con- 
crete impact on the interpretation of data on collective 
vibrations in 7-soft nuclei. In a more general perspec- 
tive, the [0(6)— U(5)]dO(5) transition of the interacting 
boson model represents a valuable theoretical laboratory 
for studying structural changes between incompatible dy- 
namical symmetries in integrable quantum systems. 
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